Asymptotic Behavior of Solutions to a Model 
System of a Radiating Gas 



Yongqin Liu* 

Faculty of Mathematics, Kyushu University 
Fukuoka 819-0395, Japan 

Shuichi KawashinW 

Faculty of Mathematics, Kyushu University 
Fukuoka 819-0395, Japan 



Abstract 

In this paper we focus on the initial value problem for a hyperbolic- 
elliptic coupled system of a radiating gas in multi-dimensional space. 
By using a time-weighted energy method, we obtain the global exis- 
tence and optimal decay estimates of solutions. Moreover, we show 
that the solution is asymptotic to the linear diffusion wave which is 
given in terms of the heat kernel. 
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1 Introduction 

In this paper we consider the initial value problem: 

I «t + E /j« + &vq = 0, ^1", t>0, 
( -Vdivg + q + Vu = 0, iGl", t > 0, 
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with the initial data 



u(x,0) = u (x), x e R n . 



(1.2) 



Here fj(u), j — 1, • • • , n, are smooth functions of u satisfying fj(u) = 0{u 2 ) 
for u — > 0, and u = u(x, t) and q — (gi, • • ■ , q n )(x, t) are unknown functions 
of x = (x\, ■ • • , x n ) G IR n (n > 1) and £ > 0. Typically, it and q represent the 
velocity and radiating heat flux of the gas, respectively. 

The system (11.11) is a simplified version of a radiating gas model in n- 
dimensional space. More precisely, in a certain physical situation, the system 
(11. ip gives a good approximation to the following system of a radiating gas, 
that is a quite general model for compressible gas dynamics where the heat 
radiative transfer phenomena are taken into account: 



where p, u, p, e and 6 are respectively the mass density, velocity, pressure, 
internal energy and absolute temperature of the gas, while q is the radiative 
heat flux, and ai and a 2 are given positive constants depending on the gas 
itself. The first three equations form the usual Euler system, which describes 
the inviscid flow of a compressible fluid, and express the conservation of 
mass, momentum and energy, respectively. We refer to the book of Courant 
and Friedrichs PQ for a detailed derivation of several models in compressible 
gas dynamics. On the other hand, the physical motivation of the fourth 
equation, which takes into account the heat radiation phenomena, is given 
in [321 15]. The simplified model (II .ip was first recovered by Hamer [8], and 
the reduction of the full system (ll.3p to (11. ip was given in [321 El EEl U\- 

There are many works on the study of the hyperbolic-elliptic coupled sys- 
tem for one-dimensional radiating gas. The earlier paper with application 
to this kind of systems is [5Tj, where Schochet and Tadmor studied the reg- 
ularized Chapman-Enskog expansion for scalar conservation laws. We refer 



to [HI [121 dH [181 HS1 E31 EH [221 [251 EH [26] for shock waves, M U EE] 



for nonlinear diffusion waves, [16] for rarefaction waves, [15] [TT| I2T] for a sin- 
gular limit and relaxation limit, and [H [291 US] for L 1 stability results, . In 
the mult i- dimensional case, Di Francesco in [2] obtained the global-in-time 
existence and uniqueness of weak entropy solutions to the system (11.11) and 
analyzed the relaxation limits. In [33J, recently Wang and Wang studied the 
initial value problem for the system (II .ip in multi-dimensions and obtained 
the pointwise estimates of classical solutions by using the method of Green 



< 



v 



p t + div(pw) = 0, 

(pu) t + div(pw <g) u + pi) = 0, 

{p(e + )} t + div{pn(e + + pu + q] = 0, 

-Vdivg + a x q + a 2 V6> 4 = 0, 



(1.3) 
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function combined with some energy estimates. Also, the stability of planar 
rarefaction waves was discussed in E] . Very recently, Ruan and Zhu [28] 
investigated the asymptotic decay rates toward the linear diffusion wave and 
also the rarefaction wave in R™, and proved the asymptotic relation of our 
Proposition 14.21 with k = and the norm L 2 instead of H s ~ 1 ~ k in the case 
2 < n < 7 by using the energy method and the semigroup argument. 

In this paper we investigate the decay rate not only to the same linear 
diffusion wave as in [28] which can be seen from Proposition I4.2| but also to 
another diffusion wave which is given in terms of the heat kernel as shown 
in Theorem 12.41 for the initial value problem ( II. ip . (I1.2p in R n for n > 2 
by applying the time-weighted energy method together with the semigroup 
argument , which removes the restriction n < 7 assumed in [23] and also 
improves their results. To this end, we first transform the system (11.11) into 
the following equivalent decoupled system (II. 4p which makes the derivation 
of our energy estimates easier, but is not essential for obtaining our main 
results, 



Then, by applying the time- weighted energy method to the decoupled system 
( 11.41) , we derive the optimal decay estimates of solutions for alln > 1 . Finally, 
using the semigroup argument, we show that the solution is asymptotic to the 
linear diffusion wave as t — > +oo, provided that n > 2. Our linear diffusion 
wave is given explicitly in terms of the heat kernel. 

The contents of the paper are as follows. In Section 2 we give full state- 
ments of our main theorems. Section 3 gives the proof of the results on the 
global existence and decay estimates of solutions. The last section gives the 
proof of the theorem on the asymptotic convergence to the linear diffusion 
wave. 

Before closing this section, we give some notations to be used below. Let 
J-~[f] denote the Fourier transform of / defined by 



and we denote its inverse transform by J 7-1 . 

L p = L p (M n ) (1 < p < +oo) is the usual Lebesgue space with the norm 
|| ■ \\lp- For 7 el, let = L*(IR n ) denote the weighted L 1 space with the 
norm 









x e R", t > o. 



(1.4) 
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Let s be a nonnegative integer. Then H s = H s (M. n ) denotes the Sobolev 
space of L 2 functions, equipped with the norm 

ii/ik- :=(£ii#/i&)*- 

k=0 

Here, for a nonnegative integer k, d k denotes the totality of all the fc-th order 
derivatives with respect to x G W 1 . Also, C k (I; H s (M. n )) denotes the space 
of /c-times continuously different iable functions on the interval / with values 
in the Sobolev space H s = H s (W n ). 

Finally, in this paper, we denote every positive constant by the same 
symbol C or c without confusion. [ • ] is the Gauss' symbol. 

2 Main theorems 

Our first theorem is on the global existence and uniform energy estimate of 
solutions to the problem (11.11) . (11.21) . 

Theorem 2.1. Let n > 1 and let s > [f] + 2 be an integer. Assume thatu G 
H s {R n ) and put E := H^oH/p- Then there exists a small positive constant 
5q such that if Eq < 5q, then the problem (11.11) . (11.21) has a unique global 
solution (u,q)(x,t) with 

u G C([0, +oo); H s (R n )), d x u G L 2 ([0, +oo); H'~\R n )), 
q G C([0, +oo); H S+1 (W 1 )) n L 2 ([0, +oo); H s+1 (« n )). 

The solution verifies the uniform energy estimate 

\Ht)\\ 2 H* + Il9(*)ll^i + f \\dMr)\\ 2 H s-i + h(T)\\ 2 Hs+ idT < CE 2 . (2.1) 

When n > 2, the solution obtained in Theorem 12. II verifies the following 
decay estimates. 

Corollary 2.2. Let n > 2 and let s > [5] + 2 be an integer. If E Q = ||mo||h s 
is suitably small, then the global solution obtained in Theorem \2. 1\ verifies the 
decay estimates 

\\d k x u(t)\\ H s- k <CE G (l+t)-i (2.2) 
for k with < k < s and 

\\%q(t)\\ H . + i-k < CE (1 + t)- k ~¥ (2.3) 
for k with < k < s — 1 . 
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The next theorem is on the optimal decay estimates of solutions for initial 
data in H s (R n ) fl L^R"). 

Theorem 2.3. Let n > 1, and let s > 3 for n = 1 and s > [f ] + 2 for n > 2. 
Assume that u G H s (W. n ) fl L 1 (M n ) and put E x : = ||^o||// s + H^olU 1 - Then 
there exists a small positive constant 5\ such that if E\ < 5\, then the global 
solution obtained in Theorem \2.1\ satisfies the decay estimates 

\\d h x u{t)\\ H s- k <C^(1 + (2.4) 

for k with < k < s — 1 and 

\\d£q{t)\\ H . + i- k < CE 1 (1 + t)-?"^ (2.5) 

for k with < k < s — 2. 

Remark 1. To obtain the optimal decay estimates of solutions stated in 
Corollary 12.21 and Theorem 12. 3[ we will use the following decay estimate for 
L°° norm of the derivative d x u as shown by N(T) in (13. 12[) : 

\\d x u(t)h°° <C(l + t)-\ 

Our final result is concerning the asymptotic profile of the global solution 
obtained in Theorem 12.31 for n > 2. First we show that for n > 2, the 
solution to the problem ( II .ip . ( II. 2p can be approximated by the solution to 
the corresponding linear problem, 

u t — Au t — Au = 0, u(x. 0) = Ur,(x), 

(2.6) 

q = -(1 - A)- x Vw. 

Then we prove that the solution to this linear problem can be further ap- 
proximated by the following simpler problem (12. 7p based on the linear heat 
equation u t — Au = 0, 

u t — Au = 0, u(x, 0) = ^0(^)1 
q = -(1 - A)~ l Vu. 
Since the solution to the linear heat equation is asymptotic to the heat kernel 

G Q (x,t) =J^- 1 [ e - | « |2< ](x) = (47rf)-te-^, (2.8) 

we thus conclude that the asymptotic profile of our global solution is given 
by the following linear diffusion wave (u*, q*)(x, t): 

u*{x,t) = MG (x,t + 1), q*(x,t) = -Vu*(x,t), (2.9) 

where M = j R „uo(x)dx denotes the "mass". The result is precisely stated 
as follows. 
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Theorem 2.4. Letn>2 and s > [f]+2. Assume thatu G # s (M n )nLj(M n ) 
and put E 2 : = H^oll^ + H m o||l}- TTten t/ie global solution (u,q) to the problem 
(11. ip . (11.21) . which is constructed in Theorem \2.3[ is asymptotic to the linear 
diffusion wave (u*,q*) in ( 12. 9 p as t ^ +oo: 

<(7J^p(0(l + *r J " 4 * 1 (2-10) 

/or < A; < s — 1 and 

11*9 - 9*)(*)lk-* < C^p(t)(l + t)-*- 4 * 3 (2.11) 
/or < k < s — 2, where p(t) = ln(l + t) for n = 2 and p(t) = 1 for n > 3. 

3 Global existence and decay estimates 
3.1 Global existence of solutions 

This subsection is devoted to the proof of the global existence result stated in 
Theorem 12.11 Since a local existence result can be obtained by the standard 
method based on the successive approximation sequence, we omit its details 
and only derive the desired a priori estimates of solutions. 

First we give a lemma which will be used in the derivation of our energy 
estimates. 

Lemma 3.1. Let 1 < p, q, r < +oo and - = - + 1. Then the following 
estimates hold: 

\\d k x (uv)\\ LP < C(\\u\\ Lq \\d k x v\\ L r + \\v\\ Lq \\d k x u\\ L r), k > 0. 

\\[d k ,u]d x v\\ LP < C(\\d x u\\ Lq \\d x v\\ L r + \\d x v\\ Lq ||<^w|| L r), k > 1. 

The proof of this lemma can be found in [20J. 

The next lemma shows the equivalence of the system (11.11) and (ll.4p . 

Lemma 3.2. The system ( 11.11) is equivalent to the decoupled system (11.41) . 

Proof. First we show that (11.11) implies (II. 4p . It follows from (ll.ll) 2 that q = 
V0 with cf) = divq — u. Therefore we see that Vdivg = VA0 = AV0 = Aq. 
Consequently, (ll.ip 2 becomes (1 — A)q + Vu = 0, which gives (11.4p 2 . Also, 
applying div to ( ll.4p 2 . we have 

divq= -(1- A) _1 A«. (3.1) 
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We substitute (13 .ip into fll.ip i and then apply 1 — A to the resulting equation. 
This yields fO]) r 

Next we derive fTTTjl from ffL4l) . We have flO) from fO]L . Apply (1 - 
A) -1 to (ll.4p i and substitute (13.11) to the result. This yields fll.ip ^ Also, we 
rewrite (13. ip as divq = u — (1 — A) -1 w and apply V, obtaining 

Vdivq = Vu - (1 - A) _1 Vw = Vw + g, 

where we have used fll.4H 2 . This gives (ll.ll) 2 . Thus the proof of Lemma [3.21 
is complete. □ 



In view of Lemma 13.21 we only need to consider the decoupled system 
(II. 4p . which makes the derivation of our energy estimates easier. As a priori 
estimates of solutions, we show the uniform energy estimate for u by using 
the equation (ll.4p i and then a similar estimate for q by making use of the 
relation (ll.4p 2 . To this end, we consider solutions u to the problem ( ll.4p i . 
(11.21) . which are defined on the time interval [0, T] for T > and satisfy 

K*)IU-<C, \\d x u(t)h°° < C (3.2) 

for < t < T. Notice that these estimates hold true if sup 0<t<T \\u{t) < C 
for s > [f ] + 2. 

Now we multiply (ll.4p i by u. A direct computation gives 

(u 2 + \Vu\ 2 ) t + 2\Vu\ 2 - 2V • {uV{u t + u + J2 fj( u )^)} 

j 

+ ^{2^) + /; wiv«i% = - E //(«)^ iv«i 2 , 

3 3 

where gj(u) = J " fj(r])r]dr]. Integrating this equality with respect to x, we 
have 

j t \\u(t)\\ 2 Hl + 2\\d x u(t)\\li < C\\d x u(t)\\L°°\\dMt)\\h- (3-3) 
To get similar estimates for the derivatives, we apply d l x to (13.11) ^ obtaining 

d l x u t - Ad x u t - Ad l x u + J2 f'j{u)d l x u Xj - V ■ J2 f^ny<i"r, 

3 3 

= -h l + V -H l , 

where U = J2j[ d L f'j( u )] u xj and H l = Y,jW x d X] , /j(u)]Vu. We multiply this 
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equation by d l x u and compute directly to get 

(\d l x u\ 2 + \Vd l x u\ 2 ) t + 2\Vd l x u\ 2 - 2V ■ {d l x uVd l x (u t + u) 

+ din £ f 3 (u)Vd l x u Xj } + ^{/;(«)(|^| 2 + \W x u\ 2 )} Xj 
j j 

= //(«K(l£«| 2 + \^ d x u ?) ~ ^h l d l x u + H l ■ Vd l x u) + 2V • (H l d l x u). 

j 

Integrating this equality with respect x and estimating the right hand side 
by applying Lemma [3. 11 we obtain 

^11^^(^)11^ -h 2||a4V^(t)||i 2 < CH^^C^IU^Iia^C*)!!^, (3.4) 

where 1 < I < s — 1. Here we have used the fact that fj(u) = 0(u 2 ),j = 
1, • • ■ , n for u — )• 0. We add (13.31) and (13 .4p for 1 < / < s — 1 and integrate 
over (0,t). This yields 



IK*)ll!r. + / \\dMr)\\ 2 H s-idr 
Jo 

< CE 2 + C [ \\d a u(T)\\ L ~\\d x u(T)\\ 2 H .-idT, 

Jo 



(3.5) 



where E = \\uq\\h s - Let 5 be a positive number (independnt of T) and 
assume that sup 0<i<T |m£)||# s < 8, where s > [~] + 2. Then the second 
term on the right hand side of (13.51) is estimated by C5 L WdxU^W^^dr . 
Therefore, choosing 5 so small that C5 < |, we arrive at the uniform energy 
estimate t 

\Ht)\\%.+ [ \\d x u(r)\\ 2 HS -idT<CE 2 . (3.6) 
Jo 

On the other hand, it follows from (I1.4IL that ||g|| HS +i < ||<9 x .w||#s-i, which 
combined with (13.61) yields 

h(t)\\W + f \\q(r)\\ 2 Ha+1 dr < CE 2 . (3.7) 
Jo 

These observations are summarized as follows. 

Proposition 3.3. Let n > 1 and let s > [|] + 2 be an integer. Assume that 
u G H s (W a ) and put E := ||wo||h s - Let (u,q) be a solution to the problem 
(11.11) . (11.21) on the time interval [0,T]. Then there is a small positive constant 
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5 independent of T such that if sup 0<t<T \\u(t) \\ H s < 5, then the solution 
verifies the following uniform energy estimate for t G [0, T]: 

\Ht)\\H* + h(t)\\W + / \\a x u(T)\\l.- 1 + \\q(r)\\ 2 Ha+l dr < CEl (3.8) 

Jo 

By virtue of the a priori estimate (13.81) for small solutions stated in Propo- 
sition 13. 3^ we can apply the continuity argument and obtain a unique global 
solution to the problem (11. ip . (II. 2p . provided that E is suitably small, say, 
Eq < Sq. The solution obtained verifies (13. 8p for t > 0. This proves Theorem 



3.2 Optimal decay estimates 

In this subsection, we obtain the optimal decay estimates of the solution 
constructed in Theorem 12.11 by using the time- weighted energy method. To 
this end, we define two time-weighted energy norms E(T) and M(T). Also, 
we introduce D(T) as the dissipation norm corresponding to E(T). 

s 

E(Tf:=J2 ^p (l + ty\\diu(t)\\l^, (3.9) 

. n 0<t<T 
J= 

D ( T ) 2 ■= E r^ + ^W^MrWns^dr, (3.10) 

„-n JO 



3=0 

s-1 



M{Tf := sup (1 + t)^\\diu(t)f H3 . 3 , (3.11) 

. n 0<t<T 
3=0 

where s > 1. To derive estimates for E(T), D(T) and M(T), we make use 
of the following time- weighted norm N(T): 

N(T):= sup (l + t)\\d x u(t)h°°. (3.12) 

0<t<T 

As for the energy E(T) and D(T), we have the following estimate. 

Proposition 3.4. Let n > 1 and s > [S] +2. Assume that u 6 if s (R n ) and 
i?o := ||^o||ff s - Then the solution to the problem (II. ip . ( II. 2p constructed 
in Theorem \2.1\ satisfies following energy estimate: 

E(T) 2 + D{Tf < CEl + CN(T)D{Tf. 
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Proof. In order to prove this proposition, it is enough to show the following 
estimates for any t G [0, T] and < j < s: 

(1 + ty\\diu(t)\\ 2 H s- 3 + Jj(l + ry\\d^u{r)\\ 2 H s-^dr 

<CEl + CN{T)D{T) 2 . 

We know from ([21]) that fl37L3]) holds true for j = 0. Now, let < k < s - 1 
and suppose that (I3.13P holds true for j = k. Then we show (I3.13P for 
j — k + 1. Multiplying (13.41) by (1 + t) k+1 , integrating with respect to t over 
(0, t) and adding for / with k + l<l<s — 1, we have 

(i + t^apMmi^-i + J?(i + r) fc+1 ||arMr)ll^- 2 ^ 

< CE 2 + CJJ(1 + r) fc ||ar^(r)||^_ fe _ lC /r 

+C/J(l + r) fe +^aXr)|Uoo||^+Vr)||^_ fe _ 1 dr. 

The second term on the right hand side is estimated by the induction hypoth- 
esis (13.131) with j = k, while the last term can be estimated by CN(T)D(T) 2 . 
Consequently, we have 

< CEl + CN{T)D{T) 2 , 

which shows that (13.131) holds true for j = k + 1. Thus, by induction, we 
have proved Proposition 13.41 □ 

By employing the optimal decay results expressed in E(T) and M(T), we 
obtain the following estimates for N(T). 

Lemma 3.5. (i) Ifn>2 and s > [§] + 2, then N(T) < CE{T). 
(ii) Ifn = 1 and s > 3, then N(T) < CM(T). 

Proof, (i) Let so = [| ]+l and 9 — By applying the Gagliardo-Nirenberg 
inequality, we see that 

<CE(T)(l + t)-U^)-^ 
= CE(T)(l + t)-fi+&. 



Since n > 2, it yields that N(T) < CE{T). 

(ii) By using the one- dimensional Gagliardo-Nirenberg inequality, we have 

<c||ax^)llill^)lli <cM{T){i + ty\ 

which gives N(T) < CM(T). This completes the proof. □ 
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Proof of Corollary [MM By virtue of Lemma 1X51 (i). we have N(T) < CE(T) 
for n > 2, which together with Proposition 13.41 gives 

E{Tf + D{Tf < CEl + CE{T)D{T) 2 . 

Put X(T) := E(T) + D(T). Then we have X(Tf < CEl + CX(Tf. This 
inequality is solved as X(T) < CEq, provided that Eq is suitably small. In 
particular, we have E(T) < CEq, which shows the desired decay estimate 
(12.21) for < k < s. Moreover, by virtue of fll.4p 2 . we have 

\\d k x q(t)\\ HS+1 - k < \\dl +l u(t)\\ HS -,- k <CE,(l + ty^, 

for < k < s — 1, where we have used (12.21) with k replaced by k + 1. This 
completes the proof of Corollary 12.21 □ 

To estimate the energy M(T), we need the following L 1 estimate of the 
solution. 

Lemma 3.6. Under the same assumptions as in Theorem \2.3\. the solution 
to the problem (II. ip . (II. 2p satisfies the following L l estimate for u: 

\Ht)\\jjL < \\u \\ L i. (3.14) 

Proof. Applying (1 — A) -1 to (ll.4l) 1 . we have 

n 

We denote by K(x) the the fundamental solution to the operator I — A, that 
is, (I-A)~ 1 u = K*u. We know that K(x) > 0, K 6 L 1 and J Rn K(x)alx = 1. 
See [2] for the details. Let js be the Friedrichs mollifier and put 

<t>s{u) = js * sign(u), ® s (u) = / <!>s(£) d t 

Jo 

We multiply (13.151) by 4>${u) and integrate the resulting equation over IR n x 
(0,t). Then, letting S — > 0, we obtain the desired L 1 estimate (I3.14p just in 
the same way as in [21 [T6l [28] . The details are omitted. □ 

By employing the time-weighted energy method together with the L 1 
estimate ( I3.14p . we get the following estimate for M{T). 
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Proposition 3.7. Let n > 1 and s > [f ] + 2. Assume t/iat u G # s (IR n ) n 
L 1 (R Tl ), and ptrf .E := ||%||# 3 an d E± := ||mo||h s + H^olU 1 - Then, if E is 
suitably small, then the solution to the problem (11. ip . (II. 2p constructed in 
Theorem \2.1\ satisfies following estimate: 

M(T) 2 < CE\{\ + JV(T))* -1 . 

Proof. In order to prove this proposition, it is enough to show the following 
estimate for any t G [0, T] and < j < s — 1: 

\\diu(t)\\l^<CEf(l + N(T)r\l + tr^. (3.16) 

For this purpose, it is sufficient to prove 

(i + ty+"\\diu(t)\\i^ + + T y+^d»^T)\\%._ j _ 1 dT 

<CEf(l + N{T)y(l + t) a -^ 

for t G [0, T] and < j < s — 1, where a > |. 

First we show (13TT7j) for j = 0. We add and ([33]) for 1 < I < s - 1, 
multiply the resulting inequality by (1 + t) a , and then integrate over (0,t). 
This yields 

(1 + tr\\u(t)r Hs + £(1 + rr\\d x u{r)\\ 2 H s-,dr 

< CEl + C/ *(l + t) - 1 ^^)!!^^ (3.18) 

+C/ '(1 + r)°||9 !e «(T)IU-IIMT)llfr.-^T =: C£ 2 + A + / 2 . 

Since < CEo by (12.11) . we can estimate the term I 2 as 

h < CE^l{l + rY\\d x u{T)f Hs ^dr < ±/ '(l + TY\\d x u{r)f Hs ^dr, 

where E is assumed to be small as CE < |. On the other hand, we divide 
ii into two parts: 

h = CJ*{1 + Ty- l \\u(r)\\l 2 dT + C/*(l + t^-IcU^II^t =: / n + J 12 . 

To estimate Ii 2 , we choose T\ so large that C(l + T X )~ 1 < |. Then we divide 
the time interval [0,t] into two parts [0,71] and [Ti,£]; here we treat the case 
t >Ti because the case t < T\ is easier. Thus we have 



I 12 < C max {(1 + tr- 1 }^ 1 \\d x u(r)\\ 2 H ^dr + \ + r)iM*OII|.-idT 
i v l t 



<CE% + J/ , (l + r)«||MT)H2r.-idr 
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where we have used ( 12. ip . Finally, we estimate the term J n . By using the 
Gagliardo-Nirenberg inequality \\u\\ L 2 < CWdxuW^WuW 1 ^ 9 with 9 = and 
applying the Young inequality, we can estimate In as 

In < Cj*(l+rr-^dXr)&u(r)\\f- d) ^ 

< l lo (! + rr\\dMr)\\hdr + CjJ(l + t)*-^\\u{t) f Ll dr 

< \ £(1 + rr\\d x u{r)\\ 2 H s-,dT + C||«o|lix(l + t) a ^, 

where we have used the L 1 estimate (13.141) and the condition a > |. Conse- 
quently, under the smallness assumption on E , we arrive at the estimate 

(1 + t) a \\u(t)\\ 2 Hs + f\l + r) a \\d x u{r)f Hs -,dr < CE%(1 + 
Jo 

which proves (I3.17P for j = 0. 

Now, let < k < s — 2 and suppose that (I3.17P holds true for j = k. Then 
we show ( I3.17P for j = k + 1. Multiplying (13 .4p by (1 + integrating 
with respect to t over (0, i) and adding up for / with k + l<l<s — 1, we 
have 

(1 + i^+lie^)!!^-! + £(! + r) fc+1+Q H^(r)||^-,- 2 rfr 

< CEl + CjJ(l + rJ^H^VrJU^dr 

+C'£(l + r) fc + 1 +"||a^(r)|| i o ||9rMr)||^- fc - 1 ^. 

Here, using the induction hypothesis (I3.17P with j = k, we can estimate the 
second term on the right hand side by CEf(l + N(T)) h (l + t) a ~ 2 . Similarly, 
we estimate the last term as 

CjJCl + rJ^IIMrJIIx-II^WII^dT 

< CN(T) jJCl+Tj^H^+Vr)!!^*- 

< CE%N(T)(1 + N{T)) k {l + t) a ~%. 
Thus we obtain 

< C£?(l + iV(T)) fe+1 (l + . 

This shows that (13.171) holds true for j = k + 1. Thus, by induction, we have 
proved Proposition 13.71 □ 
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Proof of Theorem \2.3[ By virtue of Lemma [3. 5 \ we have N(T) < C(E(T) + 
M(T)). Therefore we have from Propositions 13.41 and 13.71 that 

E(T) 2 + D{Tf < CEl + C{E{T) + M{T))D{T) 2 , 
M(T) 2 < CE\{\ + E(T) + M(T)) s ~ l . 

Put Y(T) := E(T) + D(T) + M(T). Then we have Y(Tf < CE 2 (1 + 
Y{T)) s ~ l + CY(Tf. This inequality is solved as Y(T) < CE U provided that 
Ex is suitably small, say, E\ < 5\. In particular, we have M(T) < CEi, 
which proves the decay estimate (12 .4p for < k < s — 1. Moreover, using 
(I1.4L and (12.41) with k replaced by k + 1, we obtain 

\\d k x q(t)\\ Hs+l - h < W^um^i-u < CEkil + t)-***? 
for < k < s — 2. This completes the proof of Theorem 12.31 □ 

4 Asymptotic profile 

The aim of this section is to prove Theorem 12.41 on the asymptotic profile. 
To this end, we first consider the corresponding linear problem (12. 6p . The 

l€l 2 t 

fundamental solution to (I2.6I) 1 is given by G(x,t) = J-\e ^1(4 and 
the solution to (I2.6p i can be expressed in terms of the fundamental solution 

as 

u(x,t) = (G(t) *uq)(x). (4.1) 
The solution operator G(t)* verifies the following decay property: 

Lemma 4.1. Let n > 1 and s > 0. If (j) G H s (R n ) n L 1 (R n ) ; then we have 
the following decay estimate: 

\\%G(t) * 0|| L2 < C(l + t)-S-!||0|| Ll + Ce-^\d k x 4>\\ L 2, < k < s. 

Proof By direct calculation, we have 

\\%G(t) * twit <^(4<i+4>i)iei 2fc e-^i0(oNe 

< Cf mi \tf k e-M 2t \m 2 dZ + C Jj c|>1 |e| 2fe e-'|0(O| 2 ^ 

< cj mi \e k e-^^e^\m\ 2 dc +cj m>1 \e k e-t\m\ 2 d£ 

<C(l + t)-f- fc ||0||i 1 + Ce-*||^0|li 2 . 
This completes the proof. □ 



14 



By Duhamel principle, we can express the solution of f ll.4p i (or (I3.15P ). 
fll.2p as follows: 

u(x, t) = G(t) * w - / G(t - r) * div f(u)(r)dT. (4.2) 
Jo 

Here and in the following, we use the abbreviation f{u) = • ■ ■ , f n {u)). 

We decompose the solution formula (14. 2p in the form u{t) = u(t) — F(u)(t), 
where u is the linear solution given in (14. ip and 

F(u)(x,t)= [ G(t-r) * div f{u){r)dT. (4.3) 
Jo 

First we prove that the solution to the problem (ll.4p i . (II .2p can be ap- 
proximated by the solution to the corresponding linear problem (I2.6p i . 

Proposition 4.2. Let n > 2 and s> [§] + 2. Assume that u E H s {R n ) n 
L 1 (R") and put Ei := ||«o||# 3 + H^oIIl 1 - Let (u,q) be the global solution 
to the problem (11.11) . (11.21) which is obtained in Theorem \2.1\ and let u be 
the solution to the corresponding linear problem (12.61) . which is given by the 
formula (14.11) . Then we have 

\\d k x {u - «)(t)||^-x- fc < CElp{t){\ + t)-t"^ 

for k with < k < s — 1, where p(t) = ln(l + t) for n = 2 and p(t) = 1 for 
n>3. 

Proof. Let k and m be nonnegative integers. We apply d k+m to F(u) in (14. 3 p 
and take the L? norm, obtaining 

\\d k x +m F(u)(t)\\ L2 < (/ l +/|)ll^ +m+1 G'(t-r)*/( U )(r)|| L2 dr 

By applying Lemma |4~1~| using Lemma [37X1 and noticing that fj{u) = 0(u 2 ),j = 
1, • • ■ , n for u — > 0, we have 

h <Cf*(l+t-T)-i-^\\f(<u)\\»dT 

+C^e-^\\d^f{u)\\ L ,dr 

KCj^l + t-ryi-^WuWldr 

+CQ \\u\\ L ~ \\d k x +m+l u\\ L 2dT. 
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Similarly, we have 

h < C/ t '(l+t-T)-i-5||^+ m /( M )|| L idT 

2 

+Cj t t e- t -^\\d k x + m+1 f(u)\\ L 2dr 

2 

< C/«(l + t-r)-f-2|| w || L 2||af+ m M|| L 2(ir 
+C/± e-^||w|| L oo||^ +m+1 ix|| L2 dr. 



We estimate the terms ii and ^2- Let s = [§ ] + 1 and = By using 
the Gagliardo-Nirenberg inequality and (I2.2p . we see that 



(4.5) 



\Ht)\\L- <^ll^)lli2 e ll5>(^n L 2 

< CE^l + 9 = CE 1 (1 + t)"T . 
By using ([2~2l . (g3D and (143)) . we estimate Jj as 

/i <CE%fi(l + t-T)-*- !!i ¥ i (l + T)-UT 

+CE\ $ e- 4 ?(l + r)-t-^dr (4-6) 

for fc satisfying A; + m + 1 < s, where p(t) is given above. Similarly, we can 
estimate J 2 as 

h <C7^/|(l + t-r)-T-3(l + r)-t-^dr 

+C£? j| e"^(l + r)-?-^dr (4.7) 

< C£ 1 2 p(t)(i + *)^- i ^ 

for satisfying k + m + 1 < s. 

We substitute (14. 6 p and ( 14. 7p into ( 14. 4 j) and take the sum for m with 
< ra < s — k — 1. This yields 

||fl£F(«)(t)lk-i-» < c^ 2 p(t)(i + t)-t-^ 

for with < k < s — 1. This completes the proof of Proposition 14.21 □ 

Next we consider the simpler problem for the linear heat equation ( 12.71) . 
Notice that the solution to (j2.7j) 1 is expressed as u(x, t) = (G (t) * u )(x), 
where Go(x,t) is the heat kernel given in (12.81) . By direct computation we 
see that the solution u(x, t) of (12.61) 1 is well approximated by the solution 
u(x, t) of (I2.7p i as t — > +oo: 



16 



Lemma 4.3. Let n > 1 and s>0. Assume that u G H s (R n ) n L x (]R n ) and 
put Ei = \\u \\hs + ||mo||li- Let u(x,t) and u(x,t) be solutions to (I2.6p i and 
(12.71^ . respectively. Then we have the following asymptotic relation: 

\\d k {u - u)(t)\\ H s- k < CE 1 {1 + t)-?-*- 1 , < k < s. 

Proof. Since \(G(£,t) - G (£,t)\ < C\^e~^ 2t for |£| < 1, and |(G(£,t) - 
Go(£,t)\ < Ce~ ct for |£| > 1, by direct calculation we have the following 
decay property for (G - G )(t)*: V0 G H s R L 1 , 

\\d k x (G - G )(t) * <j>\\ L * < C(l + tri-t^WWv + Ce- ct \\d k x <P\\ L2 , 
here k < s. Replace k with k + m, with iio, then we have 

R fc+m (G - G )(t) * uolb < C(l + ^-t-^-l^olUi + Ce- c '||^+ m n || L2 

< C(l + £)~*™ 5 ~ 1 ||mo||.z, 1 + Ce~ ct \\uo\\ H k+ m . 



Take sum with < m < s — k, then we complete the proof of Lemma fl~3l □ 

Finally, we show that the solution to (12. 7p is well approximated by the 
linear diffusion wave defined in (12. 9p . For this purpose, we recall the following 
well known result for the heat kernel. 

Lemma 4.4. Let n > 1 and 1 < q < 2. If (ft & L 9 (IR n ) ; then we have 

\\d k x G (t) * 0|| L2 < ertd-5)-f k > o. (4.8) 

yl/so, if (fi E L\(M. n ) and j R „ <f)(x)dx = 0, then we have 

\\d k x G {t) * <j>\\ L * <Ct~ ^U\\ L }, k>0. (4.9) 



Proof of Theorem \2.4\ We recall the linear diffusion wave (u* ,q*)(x,t) in 
(12. 9p . Since G (x, t + 1) is a solution of the linear heat equation u t — Au = 

with the initial data Go(x, 1) = (47r)~2 e~~ =: <fio(x), we see that Go(x, t + 
1) = (Go(t) * <p )(x). Consequently, we can write 

(u-u*)(t) =G {t)*{u -M<j) ), 

where M = J Rn uo(x)dx. Here 4>o(x) is a rapidly decreasing function satisfy- 
ing J Rn (f> (x)dx = 1, so that we have J Rn ( u o — M(j) Q )(x)dx = 0. Therefore, 
applying Lemma I4.4[ we deduce that 

\\d k (u - u*)(t)\\ H .-> = \\d k G (t) * (u - M<h)\\ H .-* 

<CE 2 (i + ty%-^ (4 ' 10) 
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for < k < s. Now we write u — u* = (u — u) + (u — u) + (u — u*) and apply 
Proposition 14.21 Lemma 14.31 and (I4.10p . This yields the desired estimate 
(T2TTT7]) for < k < s - 1. 

On the other hand, we have q* = —Vu* by (12.91) . which can be rewritten 
as q* = -(1 - A) _1 Vm* + (1 - A) _1 AW. Therefore, noting fOD . we find 
that 

q - q * = -(1 - A)- l V{u - u*) - (1 - A^AVw*. 
Consequently, we have 

\\d k x (Q ~ <Z*)(*)||ff.-* < \\d k x +1 (u - u*)(t)\\ Hs - 2 - k + \\d k x +3 u*(t)\\ Hs - 2 - k . 

Here, applying f)2.10p with k replaced by k + 1, we see that the first term on 
the right hand side is estimated by CE 2 p(t)(l + £)~* ~ for < k < s — 2, 
where p(t) is given in Theorem 12.41 Also, using the expression u* = MGo{t)* 
0o and applying (14.81) . we know that the second term can be estimated by 

n fc + 3 _______ 

C||mo||li(1 + t)~~~~ . These observations prove the estimate (12. lip . Thus 
the proof of Theorem 12.41 is complete. □ 
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